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We study synchronous phenomena of a coarse-grained power grid model, the swing equation, 
on small-world networks. We show that its steady state, which stands for the normal operation 
of the power systems, can be realized even if the phases are disordered. In addition, we 
prove the linear stability of steady state with small-different phases between the adjacent 
oscillators. On the other hand, a trigger of instantaneous power failure, which is described by 
the hysteretic transition, might disappear on an appropriate small-world network. This result 
suggests that the small-world connection would potentially prevent the massive blackouts. 
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1. Introduction 

The steady electric power supply and proper operation of power systems have been an 
important problem to support our daily life and avoid massive blackouts 1) . On the other hand, 
the fundamental understanding on the stable optimized operation of the power system has not 
been established sufficiently because the dynamics of the power system includes too compli¬ 
cated phenomena. 

From a point of view of the complex network analysis, the power grid is considered as a 
typical example of the large artificial networks, and its statistical properties of transmission 
networks have been investigated. Pagani and Aiello 2 ’ reviewed a lot of researches applying 
the complex network analysis to various real power grids. They showed that the real power 
grids consists of ~ 10 4 nodes, which connect with edges of transmission lines, and their 
mean degrees, the mean numbers of edges of transmission lines per a node, are in the range 
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of 2 ~ 4, respectively. If the mean distance is small and the cluster coefficient is large the 
network becomes a small world, which has been considered as a distinctive character of 
natural networks 3 '. However, in the case of the power grids, the mean distance is small but 
the cluster coefficient variously changes. In other words, it has not been clear whether there 
are benefits of small world for the power grid network. 

In the present study we discuss the stability and transition of power systems on the small- 
world network generated on the basis of the Watts-Strogatz model 4 '. The model creates the 
small world by the random reconnection of the regular loop. We deal with the swing equa¬ 
tion, which is a coarse-grained model on the phases of the alternating currents in the power 
system 5 ' 6 ', on the Watts-Strogatz model. It should be noted that the swing equation without 
the inertia term corresponds to the Kuramoto model 7 ', which is a solvable oscillator model 
of synchronous phenomena. The swing equation has thus been attracting many researchers’ 
interest not only as a model of power grids but also as an extended problem of Kuramoto 
oscillators. Indeed the swing equation has been analyzed by Tanaka et al. 8 ' in particular on its 
phase transition on the complete graph for various torque distributions. They found that the 
oscillators have the first-order hysteretic phase transition, the process of which depends on 
the type of the torque distributions and the mean-field approximation gives good predictions 
on the phase transition. The synchronous conditions of the swing oscillators was studied by 
Dorfler and Bullo 9 ' through the singular perturbation from the Kuramoto oscillators. They 
proved that the oscillators converge to the phase-lock solution of Kuramoto oscillators when 
the mass of swing oscillators is sufficiently small. The simplicity of the swing equation drives 
other researches such as the distributed power generating system 10 ' and the stability of steady 
states on the real power networks 11 '. 

In order to compare the swing equation to the Kuramoto oscillators we briefly review 
the studies on the Kuramoto oscillators. The phase transition of the Kuramoto oscillators on 
the Watts-Strogatz network was studied by Hong et al. 12 '. They reported that the oscillators 
appear the second-order phase transition with the frequency synchronization similar to that 
on the complete network. On the complete graph, the process of phase transition depends 
on the type of frequency distribution 7 ', which differs from the case of random graph. Um et 
al. 13 ' studied the phase transition of Kuramoto oscillators on the random graph varying the 
type of the frequency distribution. They found that since the random connections violate the 
transition characters the oscillators have the common second-order phase transition regardless 
of the frequency distributions. 

In this paper, we discuss the stability and phase transition of the swing equation on the 
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Watts-Strogatz model to study properties of the power systems on the small-world network. 
Although the analysis on the small world following Watts-Strogatz model is not enough to 
directly understand the real power grid, we focus on the stability and transition of swing 
equation on the Watts-Strogatz model as a simplified model of the real power grid. 

The remaining of the paper is organized as follows. The second section fixes the nota¬ 
tion and definition of the swing equation and the Watts-Strogatz model. In the third section 
we prove the linear stability of the steady state with the small-different phases between ad¬ 
jacent oscillators. This stability allows a different type of the steady state from that with the 
well-known synchronization. In the forth section we demonstrate several numerical compu¬ 
tations and report the following results. In the case of the random graph, the oscillators have 
the common first-order hysteretic phase transition between the ordinary steady state with 
synchronization and the unsteady state regardless of the type of the torque distribution. In 
the case of small-world network, on the other hand, the oscillators become the steady states 
with small-different phases without synchronization, and the hysteretic transition seems to 
disappear on an appropriate small-world network. This fact implies that on the small-world 
network we might be capable to operate the power system without fear of the blackout. The 
fifth section is devoted to conclusion and discussion. 

2. Problem setting 

We deal with the dimensionless swing equation on a small-world network. The phase of 
/th-gcncrator/consumcr, 0 ; - e [0,2 n), is governed by 

N 

mOi + 6i -tOj + K ^ Ajj sin (6j - 0,), (1) 

j 

where 0 ; - = d 2 djldt 2 and 0,- = ddj/dt are the time derivatives. Here t denotes the time, and 
N represents the total number of generators and consumers 56) . The steady state, 0, = 0(7 = 
1, • • ■ N), represents the normal operation in the context of the power systems. On the other 
hand the unsteady state, 0, ± 0, signifies the blackout. In addition m and K are the mass and 
the coupling constant, respectively. The torque follows a distribution g(a>) with a vanishing 
mean. Then tOj > 0 indicates the generator providing energy to the power systems, and oj, < 0 
stands for the consumer expending the energy transferred from the generators. It is noted that 
the zero mean of the torque distribution is necessary for the existence of steady states. The 
second term on the right-hand side of the equation (1) indicates the energy transfer between 
each node, generator and consumer, and A, is the (7, j) element of the adjacency matrix A. 
If the node i and j are connected A , 7 = 1 and otherwise A,-j = 0. The network is assumed to 
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be undirected with neither loop nor multiedge. Then, A is a symmetric matrix without any 
diagonal elements, A, 7 - = A /V and A„ = 0. The degree of z'th node is given by k t - and 

the mean degree is obtained by (k) = X, kj/N. We note that, from a point of view of the phase 
reduction theory, the inertial term, mO,, is ad hoc but it is necessary for describing abrupt 
phenomena on the power grids such as the instantaneous power failure 15 ^ 

In order to investigate the behavior of the phase following the swing equation on a typical 
small-world network, we employ the Watts-Strogatz model 4) . The model can generate various 
small-world networks by modifying the regular network with the reconnection probability P. 
In the case of P - 0 the network is regular, and the mean distance and cluster coefficient are 
large. On the other hand, in the case of P = 1, the network changes into the random graph with 
the short mean distance and small cluster coefficient. In the range of 0 < P < 1 the network 
becomes small world; the small mean distance and large cluster coefficient. Notice that, on the 
real power grids, the mean distance is small but the cluster coefficient can variously change 2) . 

3. Linear stability of steady states with small-different phases 

In this section, we first prove the linear stability of a steady state with small-different 
phases between adjacent oscillators on an undirected network, 

N 

0 = tOi+K Yj a U sin( e°j - 6 ?) i = 1, • • • , N, 

7=1 

\e () t - e ( )\ < ^ for {(i,j)\Aij = 1}. 

Substituting = 0° + 0, exp(crf) to (1) and neglecting the second order term of 6 n we obtain 
a linearized equation of 0°, 

N 

cr 2 mdj + crdj = K ^ A tj cos(^ - 0°)(0 7 - 0,-), 

7=1 

where " is omitted for simplicity. Taking a product with 0, over i leads to the energy equation 
of infinitesimal disturbance, 



Here we use the fact A, 7 = A /( . Since the equation (2) is an algebraic equation of the second 
degree with respect to cr, we obtain the roots 
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where V and P are the kinetic energy and potential, 

v=l Ji'L e ' >0 - p = V Z lA'costs 1 ’-el'Kfl,-fl f ) 2 > °. 

i= 1 i= 1 7=1 

Therefore Re[cr ± ] < 0. We conclude that a steady state with small-different phases between 
adjacent oscillators on an undirected network is linearly stable. It should be emphasized that 
the steady state does not necessarily imply the phase synchronization, although the steady 
state with small-different phases includes the ordinary phase lock state. Dorfler and Bullo 9) 
analytically studied the phase-lock state of the swing equation using the singular perturbation 
from the Kuramoto oscillator (m = 0), and they proved that, if m is sufficiently small and 
NK > max (m, - ojf), the oscillator asymptotically converges to a phase-lock solution of 

I <(< j<N 

Kuramoto oscillator with the order of O(m). Notice that our result is free from the limitation 


The remaining problem is where such a steady state exists. In the next section we show 
several numerical results to demonstrate that this different type of the steady state with small- 
different phases, the linear stability of which is here proven, appears as the reconnection 
probability gradually decreases. 


4. Phase transition on a small-world network 

In the present section we perform numerical integration of the swing equation and show 
the phase transition of the swing oscillators on the Watts-Strogatz network. We then take three 
types of the torque distributions; the uniform distribution, Gaussian and Bimodal Gaussian, 


8u(oj ) = 


1/2 M < 1 
0 M > 1 


= -p= exp 
Yljt 

8bg(^>) = 2^Sg(^ - 1 ) + 8g(m + !))• 

The numerical integration is performed by the second-order Adams-Bathforth method. We 
take the time step 8t = 0.05 and increase the end time of integration up to Tj = 10 4 to 
reach the steady state. The initial condition of phases and angular velocities is generated from 
uniform distributions in the range of [-1/2,1/2]. In order to evaluate the phase transition of 
the swing oscillators we estimate the order parameter, R, which characterizes the phase-lock 
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state. 


* = ir ly 


exp (idi(t)) 


dt, 


where T s = 200 and T w = T f - T s . When the oscillators are in a phase-locked state. A 1 - 1. 
On the other hand, in the disordered state, R - 0(1/N). In order to check the steadiness of 
the state we compute the squared angular velocity, 




dt. 


When the squared angular velocity vanishes, the steady state, the normal operation of the 
power systems, is realized. On the other hand the unsteady state with nonzero squared angular 
velocity signifies the blackout. In the case of K = 0, the phase states are given by 0,-(t) = 
cotf + Cn exp (-t/m) + Cq where C n and Cq are constants depending on the initial condition. 
If K is small the squared angular velocity is about the variance of the torque distribution, 
f ci) 2 g(to)dcL>. The number of oscillators is up to N = 10 4 . We take the ensemble average over 
100 samples of the randomly generated networks and torques. 

Figure 1 shows the order parameter for each torque distribution, while we plot the squared 
angular velocity in Fig. 2 in order to check steady states. We here show the representative 
results of the case m = 5 and (k) = 4. We remark that we have confirmed that the results do 
not qualitatively change with varying m and (k). In addition the results do not qualitatively 
change by varying the end time and the system size. 

In the case of the random graph, P = 1.00, as the coupling constant increases, the or¬ 
der parameter increases discontinuously and the oscillators become the phase-lock state with 
R ~ 1. On the other hand, when we gradually decrease the coupling constant, the order param¬ 
eter falls discontinuously and the disordered state appears. As the order parameter changes 
discontinuously between 0(1/N) and 0( 1), the squared angular velocity also changes discon¬ 
tinuously between f oj 2 g(a>)doj and 0(0) as in Fig. 2. The coupling constant of transition from 
the disordered unsteady state to phase-lock steady state is larger than its inverse. Thus, the os¬ 
cillators have the first-order hysteretic phase transition in the case of the random graph. Even 
if the torque distribution is different the transition process seems quantitatively unchanged. 

On the other hand, on the small-world networks, which is generated by gradually decreas¬ 
ing the reconnection probability, the transition process becomes smooth compared to the case 
on the random graph. In the region with the high coupling constants the squared angular ve¬ 
locity becomes as small as that of the random graph. In particular let us focus on the case with 
P = 0.01. The order parameter hardly varies from 0(1/N), while, as the coupling constant 
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coupling constant 



P=1.00 '—|-' P=0.I0 X ' P=0.01 —Q- 



Fig. 1. (Color online) Order parameter ( N = 10 4 , Tf = 10 4 ): (a) uniform, (b) Gaussian and (c) bimodal 
Gaussian torque distributions. The horizontal and vertical axes indicate the coupling constant and the order 
parameter/?, respectively. The red crosses, blue saltires and green circles indicate the cases of P — 1.00,0.10 

and 0.01, respectively. The error bars are the standard deviation of realizations. 
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Fig. 2. (Color online) Same as fig. 1 but the squared angular velocity ( N = 10 4 , Tf = 10 4 ). 


changes, the squared angular velocity changes between f u> 2 g(to)dco and (9(0) through the 
different paths. Notice that the squared angular velocity exhibits the smooth hysteretic transi¬ 
tion. This fact ensures that the steady state can be realized although the order parameter is a 
small value. As proven in the previous section, the steady state with small-different phases is 
linearly stable. In other words, the small-world network admits the steady state different from 
the ordinary phase-lock state. 

In addition, the transition processes of phases and squared angular velocity are clearly 
different from that of the random graph. We have to emphasize that the transition point of 
the squared angular velocity is different from that of the order parameter. Although the or¬ 
der parameter is often investigated to identify the nontrivial phase transition, the steadiness 
is also important to characterize the behavior of the oscillators. In the context of the power 
system, the steadiness should be rather highlighted since the steady state represents the nor¬ 
mal operation. In this sense, the hysteretic transition should be avoided in order to keep the 
normal operation on the small-world network. Let us closely investigate the existence of the 
hysteretic transition in the squared angular velocity below. 

We compute a rate of transition samples in order to investigate smooth hysteretic process 
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Fig. 3. (Color online) Rate of transition^ = 10 4 , Tf = 10 4 ): (a) uniform, (b) Gaussian and (c) bimodal Gaus¬ 
sian torque distributions. The top and bottom columns show increasing and decreasing processes, respectively. 
The maroon squares, red circles, green triangles, blue-reverse triangles and navy rhomboids indicate the cases of 
P = 1.00,0.50,0.10,0.05 and 0.01, respectively. The horizontal and vertical axes indicate the coupling constant 
and the ratio of transition, respectively. 


of the squared angular velocity in detail. We set the threshold which judges a sample in the 
steady state or not at f a> 2 g(aj)dio/2, which is the half of the squared angular velocity at 
K = 0. As the coupling constant increases/decreases we count the number of samples of 
which the squared angular velocity is more/less than this threshold. The rate of transition 
is given by dividing this number of counted samples by the number of all samples. Figure 
3 shows the rate of transition for each torque distribution. The rate of transition is a near 
step function, and the critical points of transition could be distinguished clearly. We find 
that as the reconnection probability decreases the difference between two critical coupling 
constants decreases until P - 0.1, and as the reconnection probability further decreases this 
difference increases monotonically. It should be noticed that on Fig. 2 in the case of P = 0.1 
the difference of the hysteretic paths is smaller than that of the cases of P - 1.0 and 0.01. This 
tendency is consistent with the results of the rate of transition. These results imply that for 
any torque distribution the hysteretic process of the squared angular velocity would change 
and possibly disappear around P - 0.1. It is noted that we have obtained this result using 
10 4 samples and N = 500. We restrict ourselves to just saying possibly that the hysteretic 
transition disappears on some small-world network. 

In addition, we study the steady state at the extremely high coupling constant and find 
that, as the reconnection probability decreases, the order parameter decays but its degrade is 
not significant (not shown). This result suggests that the phase-lock state is realized regardless 
of the reconnection probability as long as the coupling constant is sufficiently large value. We 
also study the correlation of the phases of steady oscillators, and we find that the oscillators 
distant about the mean distance of network are correlated (not shown). It is noted that, as 
the reconnection probability decreases, the mean distance increases. This result implies that 
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the correlation between distant oscillators would lead to the variety of phases of steady state. 
Especially, the steady state without the phase synchronization would arise from the slightly 
longer distant correlation rather than that of the random graph, such as the case of P = 0.01 
shown in Figs. 1 and 2. 

5. Conclusion and discussion 

We study the phase transition of swing equation on a small-world network. We prove the 
linear stability of the steady states with small-different phases between adjacent oscillators. 
This result permits the steady state without the phase lock. 

Next, we study the phase transition on the Watts-Strogatz network for various torque 
distributions. In the case of the random graph, we find that the oscillators have the lst-order 
hysteretic phase transition. The transition process seems to be quantitatively unchanged by 
the type of torque distributions. On the other hand Um et al. l3) found that the Kuramoto 
oscillators on the random graph have the common 2nd-order phase transition regardless of the 
type of the frequency distributions. This result suggests that on the random graph the uniform 
phase transition independent of the distribution is a common property between the swing and 
Kuramoto oscillators. In addition, we have employed the mean field approximation 14) but the 
approximation does not work well. Um et al. 13) reported that the mean field approximation 
for the Kuramoto oscillators on the random graph is not consist with the numerical results. 
Therefore this result implies that the earlier mean field method could not deal with the random 
distributed interactions. 

In the case of the small-world network, we find that as the reconnection probability 
decreases the transition process of the squared angular velocity becomes discontinuous to 
smooth. In spite of incoherent phases the oscillators become the steady state as the coupling 
constant increases. Moreover, at the marginal reconnection probability, the hysteretic transi¬ 
tion of the squared angular velocity seems to disappear. Hong et al. I2) studied the Kuramoto 
oscillators on the Watts-Strogatz network and found that the 2nd-order phase transition of os¬ 
cillators accompanies the frequency synchronization. These results implies that on the small- 
world network the transition of swing oscillators would be different from that of Kuramoto 
oscillators. 

The hysteretic phenomena have been considered as the triggers of the instantaneous 
power failures 511153 . In the present study the hysteresis of the swing equation seems to dis¬ 
appear provided the network has appropriate small-world connections. This result implies 
that the power grid system would be more stable if the transmission network is built up 
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with small-world connections. This expectation coincides with the suggestion that the small- 
world connections of the real networks are optimized results of the stability and basin of 
synchronization 16) . Since our coarse-grained model is homogeneous coupling and no phase 
shift term, we neglect the inhomogeneous capacities of transmission lines and reactive power. 
Those factors vary the synchronous properties of Kuramoto oscillators 17) , and therefore more 
detail investigations are required in future. 
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